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Abstract . This report is the third of a series (Refs. 1-4) 
dealing with the determination of optimal trajectories for the aero- 
assisted flight experiment (AFE) . The intent of this experiment is to 
simulate a GEO-to-LEO transfer, where GEO denotes a geosynchronous 
Earth orbit and LEO denotes a low Earth orbit. Specifically, 
the AFE spacecraft is released from the space shuttle and is 
accelerated by means of a solid rocket motor toward Earth, so as 
to achieve atmospheric entry conditions identical with those of 
a spacecraft returning from GEO. During the atmospheric pass, 
the angle of attack is kept constant, and the angle of bank is 
controlled in such a way that the following conditions are satisfied: 
(a) the atmospheric velocity depletion is such that, after 
exiting, the AFE spacecraft first ascends to a specified apogee and 
then descends to a specified perigee; and (b) the exit orbital 
plane is identical with the entry orbital plane. The final maneuver , 
not analyzed here, includes the rendezvous with and the capture 
by the space shuttle. 

In this report, the trajectories of an AFE spacecraft are 
analyzed in a 3D-space, employing the full system of 6 ODEs describing 
the atmospheric pass. The atmospheric entry conditions are given, 
and the atmospheric exit conditions are adjusted in such a way that 
requirements (a) and (b) are met, while simultaneously minimizing 
the total characteristic velocity, hence the propellant consumption 
required for orbital transfer. Two possible transfers are 
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considered: (IA) indirect ascent to a 178 NM perigee via a 197 NM 

apogee; and (DA) direct ascent to a 178 NM apogee. 

For both transfers, two cases are investigated: (i) the 

bank angle is continuously variable; and (ii) the trajectory is 
divided into segments along which the bank angle is constant. 

For case (ii) , the following subcases are studied: two segments, 

three segments, four segments, and five segments ; because the time 
duration of each segment is optimized, the above subcases involve 
four, six, eight, and ten parameters, respectively. 

A surprising result of the analysis is that the optimal 
trajectories of cases (i) and (ii) coalesce into a single trajectory: 
a two-subarc trajectory, with the bank angle constant in each 
subarc (bang-bang control) . Specifically, the bank angle is 
near 180 deg in the atmospheric entry phase (positive lift projection 
phase) and is near 0 deg in the atmospheric exit phase 
(negative lift projection phase). Another surprising result is that, 
during the atmospheric pass, the peak values of the changes of the 
orbital inclination and the longitude of the ascending node are 
nearly zero; hence, the peak value of the wedge angle (angle between 
the instantaneous orbital plane and the initial orbital plane) is 
nearly zero. This means that the motion of the spacecraft is 
nearly planar in an inertial space. 

The guidance implications of the above results are discussed. 

Key Words. Flight mechanics, hypervelocity flight, atmospheric 
flight, optimal trajectories, aeroassisted flight experiment, 
aeroassisted orbital transfer, guidance. 
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1 . Introduction 

In previous reports (Refs. 1-2), the equations of motion of 
an AFE vehicle were derived first in an Earth-fixed system (Ref. 1) 
and then in an inertial system (Ref. 2). In this report (Ref. 3), 
the problem of the optimal trajectories is formulated, solved, and 
analyzed. In a subsequent report (Ref. 4), more technical details 
are provided in the form of data, tables, and graphs. 

The aeroassisted flight experiment (AFE) refers to the study 
of the free flight of an autonomous spacecraft, shuttle-launched 
and shuttle-recovered. Its purpose is to gather atmospheric entry 
data for use in designing aeroassited orbital transfer vehicles 
(AOTV) . The intent of this experiment is to simulate a GEO-to-LEO 
transfer, where GEO denotes a geosynchronous Earth orbit and LEO 
denotes a low Earth orbit (Ref. 5) . 

In an actual AOT tranfer, GEO-to-LEO, the maneuver is initiated 
with a propulsive impulse at GEO so as to decelerate the spacecraft 
and force it into an elliptical transfer orbit leading to atmospheric 
entry. In a simulated AOT transfer, GEO-to-LEO, the maneuver is 
initiated by releasing the spacecraft from the space shuttle, 
which is flying at the altitude of 160 NM above the Earth surface. 

By means of a solid rocket motor, the AFE spacecraft is accelerated 
toward Earth, so as to achieve atmospheric entry conditions 
identical with those of a spacecraft returning from GEO. Thus, 
for the purposes of this report, the atmospheric entry conditions are 
to be considered as given. 
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During the atmospheric pass, the angle of attack of the AFE 
spacecraft is kept at the constant value a = 17 deg. This value 
is such that the lift L is negative. Control of the AFE spacecraft 
is achieved via the angle of bank y. Hence, the projected lift 
Lp = Lcosy (the vertical component of the lift vector) can be 
made positive or negative at will, depending on the value of y. 

The time history of the control, the bank angle y(t), is 
subject to two requirements: (a) the atmospheric velocity depletion 

must be such that, after exiting, the AFE spacecraft first ascends to 
a specified apogee and then descends to a specified perigee; and 
(b) the exit orbital plane must be identical with the entry orbital 
plane; this is the same as stating that the exit value of the 
wedge angle must vanish; see Ref. 5. Here, one additional 
requirement is considered: (c) the propellant consumption required 

for orbital transfer must be as small as possible, implying 
that the characteristic velocity AV must be as small as possible; 
in turn, minimizing AV is equivalent to maximizing the exit value 
of the horizontal component of the inertial velocity; 
this implies either maximizing the exit value of the inertial 
velocity or minimizing the exit value of the inertial path 
inclination. 

In this report, the optimal trajectories of the AFE spacecraft 
are determined with reference to a 3D-space and employing the full 
system of 6 ODEs describing the atmospheric pass. In the general 
formulation of Refs. 1-2, the effects due to the rotation of the 
Earth and the oblateness of the Earth are included. In this report 
and the subsequent report (Refs. 3-4), the former are included. 
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but the latter are excluded, since they have been found to be small. 
The resulting optimal control problem is solved by means of the 
sequential gradient-restoration algorithm (SGRA, Refs. 6-12). 

While this algorithm is available in both the primal formulation 
(Refs. 6-9) and the dual formulation (Refs. 10-12) , the former 
is used in this report, based on previous experience with various 
AOT problems (Refs. 13-19) . 

Using SGRA, the optimal trajectories are computed for two 
possible transfers: (IA) indirect ascent to a 178 NM perigee via 

a 197 NM apogee; and (DA) direct ascent to a 178 NM apogee. For 
each transfer, two cases are investigated: (i) the bank angle 

is continuously variable; hence, y (t) is treated as a control; 
and (ii) the trajectory is divided into segments along which the 
bank angle is constant; hence, for each segment, y is treated as 
a parameter. For case (ii) , the following subcases are studied: 

2, 3, 4, 5 segments; because the time duration of each segment is 
also a parameter to be optimized, the above subcases involve 4, 6, 

8, 10 parameters, respectively. 

For comparison purposes and only for Transfer (IA) , a reference 
trajectory is also considered: this is a 5-segment trajectory, 
close to the nominal trajectory given in Ref. 5. 

To sum up, this report considers ten optimal trajectories, 
five for Transfer (IA) and five for Transfer (DA) . These optimal 
trajectories are compared with (and are found to be superior to) ' 
the reference trajectory in terms of the main quantities of 
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interest, namely: total characteristic velocity, peak heating 
rate, peak dynamic pressure, peak change of orbital inclination, 
peak change of the longitude of the ascending node, and peak 
value of the wedge angle. 

The report ends with a discussion of the guidance implications 
of the results obtained on optimal trajectories. 

Content. Section 2 contains the notations. The system 
description is given in Section 3, and the optimization problem 
is formulated in Section 4 . The experimental data are given in 
Section 5, and the numerical results are presented in Section 6. 

The guidance implications of the results are discussed in Section 7 , 
and the conclusions are given in Section 8. 
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2 . Notations 

Throughout this report, the following notations are employed 
C D = drag coefficient; 

C T = lift coefficient; 

Xj 

D = drag, N; 

2 

g = local acceleration of gravity, m/sec ; 
h = altitude, m; 

i = orbital inclination, rad; 

L = lift, N; 
m = mass, kg; 

r = radial distance from the center of the Earth, m; 
r g = radius of the Earth, m; 

r = radius of the outer edge of the atmosphere, m; 

cl 

2 

S = reference surface area, m ; 

t = time, sec; 

V = velocity, m/sec; 

V = circular velocity at r = r , m/sec; 

a & 

a = angle of attack, rad; 
y = path inclination, rad; 

n = wedge angle, rad; 

0 = longitude, rad; 

y = bank angle, rad; 

3 2 

y e = Earth's gravitational constant, m /sec ; 

3 

p = air density, kg/m ; 
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t = final time, sec; 

4> = latitude, rad? 

X = heading angle, rad; 

w = angular velocity of the Earth, rad/sec; 

n = longitude of the ascending node, rad; 

AV = characteristic velocity, m/sec. 

Subscripts 

0 = entry into the atmosphere; 

1 = exit from the atmosphere; 

11 = perigee for Transfer (IA) or apogee for Transfer (DA) ; 

22 = apogee for Transfer (IA) . 

Superscripts 

• = derivative with respect to time; 

= quantity defined in an inertial system. 


Acronyms 


AFE 

= 

aeroassisted flight experiment; 

AFEV 

= 

aeroassisted flight experiment vehicle; 

AOT 

= 

aeroassisted orbital transfer; 

AOTV 

= 

aeroassisted orbital transfer vehicle; 

DA 

= 

direct ascent; 

DP 

= 

dynamic pressure, N/m^ ; 

GEO 

= 

geosynchronous Earth orbit; 

HEO 

= 

high Earth orbit; 

HR 

= 

, 2 

heating rate, W/m ; 

IA 

s 

indirect ascent; 
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LAN 

= 

longitude of the ascending node, rad; 

LEO 

= 

low Earth orbit; 

ODE 

= 

ordinary differential equation; 

OT 

= 

optimal trajectory; 

RT 


reference trajectory; 

SGRA 

= 

sequential gradient-restoration algorithm; 

WA 

= 

wedge angle, rad. 
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3 . System Description 

The motion of the AFE spacecraft takes place partly in the 
atmosphere and partly in space. For the purposes of this report, 
the initial point is the atmospheric entry point; it corresponds 
to GEO return conditions and is fixed; the final point is located 
at the altitude of 178 NM and corresponds to circularization into 
LEO. We consider two transfer maneuvers; (IA) indirect ascent 
to a 178 NM perigee via a 197 NM apogee; and (DA) direct ascent 
to a 178 NM apogee. 

For Transfer (IA) , the key points of the maneuver are 
these: point 0, atmospheric entry; point 1, atmospheric exit; 
point 22, apogee (h = 197 NM) of the first post-atmospheric 
transfer orbit; and point 11, perigee (h = 178 NM) of the second 
post-atmospheric transfer orbit. Propulsive impulses are applied 
at two points: at the apogee 22 in order to raise the height of the 
perigee; and at the perigee 11 in order to circularize the motion. 

For Transfer (DA), the key points of the maneuver are these: 
point 0, atmospheric entry; point 1, atmospheric exit; and point 
11, apogee (h = 178 NM) of the post-atmospheric transfer orbit. 

A propulsive impulse is applied at only one point: at the apogee 11 
in order to circularize the motion. 

For the atmospheric portion (h < h a ) of the trajectory of 
the AFE spacecraft, we employ an Earth-fixed system; for the space 
oortion of the trajectory (h > h ), we employ an inertial system; 

— a. 

here, h = 400,000 ft = 121.9 km denotes the thickness of the 

3 
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atmosphere. For h < h , we compute the air density using the us 
Standard Atmosphere, 1976 (Ref. 20); for h > h , we assume that the 

~ ^ 3 . 

air density is zero. For both the atmospheric portion and the space 
portion of the trajectory, we neglect the effects due to the 
oblateness of the Earth; we assume that the gravitational field 
is central and obeys the inverse square law. 

Atmospheric Pass . With reference to the atmospheric portion 
of the trajectory of the AFE vehicle, the following additional 
hypotheses are employed: (a) the atmospheric pass is made with 

engine shut-off; hence, the AFE spacecraft behaves as a particle 
of constant mass; (b) the angle of attack is constant, a = 17 deg; 

(c) under extreme hypersonic conditions, the dependence of the 
aerodynamic coefficients on the Mach number and the Reynolds number 
is disregarded; (d) the sideslip angle is zero; hence, the side force 
component of the aerodynamic force is zero. The only control 
is the angle of bank y. 

Differential System . With the above assumptions and upon using 
an Earth-fixed system, the equations of motion include the 
kinematical equations (Ref. 1) 

« 

0 = Vcosy cosy/rcos<{> , (la) 

<j> =-Vcosy siny/r, (lb) 

r = Vsiny, (lc) 

and the dynamical equations (Ref. 1) 

V =-D/m - gsiny 

2 2 

+ w r(siny cos <j> + cosy siny cost)) sin<|>) , 


(Id) 
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y = (L/mV) cosy + (V/r - g/V)cosy + 2u> cosy cos<J> 

2 2 

+ ( co r/V) (cosy cos <j> - siny siny costf> sin<£) , (le) 

X = (L/mV) siny/cosy + (V/r)cosy cosy tan<}> 

+ 2w(sin<f> + tany siny cos<f>) 

+ (a) r/V) cosy cosifi sintf>/cosy. (If) 

In the dynamical equations, the symbol w denotes the angular 
velocity of the Earth? terms linear in w are due to the Coriolis 
acceleration; terms quadratic in w are due to the transport 
acceleration. Also in the dynamical equations, the acceleration 
of gravity is given by 

g = H e /r 2 , (2) 

where ]i & denotes the Earth's gravitational constant. In addition, 
the aerodynamic forces are given by 

D = (l/2)C D (a)p(h)SV 2 , (3a) 

L = (l/2)C L (a)p(h)SV 2 , (3b) 

where the air density p depends on the altitude h, with 

h = r - r . (3c) 

e 

Since a is constant, both the drag coefficient and the lift 
coefficient are constant. The control of the spacecraft is 
the angle of bank y . Equations (1) must be integrated, subject 
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to (2) and (3), over the time interval 0 <_ t £ x . Here, the 
initial time t = 0 corresponds to atmospheric entry; and the 
final time t = t corresponds to atmospheric exit. 

Transformation Relations . The following transformation relations 
allow one to pass from quantities computed in an Earth-fixed 
system to quantities computed in an inertial system, and viceversa: 


0 = 0+ cot, (4a) 
J = <j>, (4b) 
r = r, (4c) 
Vcosy cosy = Vcosy cosy + ur cos<f>, (4d) 
Vcosy siny = Vcosy siny, (4e) 
Vsiny = Vsiny. (4f) 


The first three equations are linear and refer to the state 
variables appearing in the kinematical equations. The 
next three equations are nonlinear and refer to the state variables 
appearing in the dynamical equations; in spite of the 
nonlinearity, these equations can be solved explicitly to yield 
V, y, y in terms of V, y, y, and viceversa. This step is omitted 
for the sake of brevity? see Refs. 1-2. 
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Orbital Elements . Once the state variables are known in 
the inertial system, one can compute some important orbital 
quantities such as the orbital inclination i and the longitude of 
the ascending node ft. These quantities are supplied by the 
relations 


cosi = cos(f> cosy/ 

(5a) 

sin(0 - ft) = coti tan<}> . 

(5b) 


Another important quantity is the wedge angle g, which is the 
angle between the instantaneous orbital plane and the entry orbital 
plane. This angle is given by 

cosn = sini sini^ cos (ft - ftp) + cosi cosi^. (5c) 

Summary . To sum up, in the Earth-fixed system, the equations 
governing the atmospheric pass include the differential system 
(1) and the analytical and functional relations (2) -(3). In 
this formulation, the independent variable is the time t, 0 <_ t <_ x . 
The dependent variables include six state variables [0(t), 4>(t), 
r(t), V(t), y (t) , x(t)]f one control variable [y(t)], and one 
parameter (t) . If the initial values of the state variables are 
given and if the bank angle program is prescribed, the system 
can be integrated in forward time. Then, using (4) , one can 
convert the state variables computed in the Earth-fixed system 
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into state variables computed in the inertial system. Finally, 
using (5) , one can determine the orbital inclination, the 
longitude of the ascending node, and the wedge angle. 

Initial Conditions . At atmospheric entry, the initial values 
of the state variables 0 Q , 4> Q , r Q = r & , V Q , y Q , x 0 are given in 
the inertial system. Hence, the initial values of the orbital 
elements 1 q, Qq are known. By definition, the initial value of 
the wedge angle satisfies the relation = 0. Using the 
transformation relations (4), the initial values of the state 
variables e Q , <t> Q , r Q = r a , V Q , y Q , x 0 become known in the Earth- 
fixed system. 

Final Conditions . The final time x is free and is to be 
determined indirectly as the time instant at which the spacecraft 
exits from the atmosphere; hence. 


r = r . (6a) 

1 a 

Application of energy conservation and angular momentum 
conservation to the exit— to— apogee transfer orbit yields either 
the relation 


(IA) 


(2V^ - V*) - 2r 22 r a v; + rjv'cos*^ = 0 


~2 


2 ~ 2 


22 


or the relation 


(DA) 


ril^a 




~2 2~2 2 ~ 

2r . , r - v = + r V cos y 
11 a a a 1 1 


(6b) 


0. 


(6c) 
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Equation (6b) applies to the indirect ascent case; Equation (6c) 
applies to the direct ascent case. These relations guarantee 
that, after exiting, the spacecraft ascends to the specified apogee. 
Finally, the orbital plane at atmospheric exit must be identical 
with the orbital plane at atmospheric entry; hence, the exit value 
of the wedge angle must vanish, 

= 0. (6d) 

Use of the transformation relations (4) in conjunction with (5) 
allows one to transform the final conditions (6) to include 
quantities computed in the Earth-fixed system. 
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4 . Optimal Control Problem 

Subject to the previous constraints, different optimal control 
problems can be formulated, depending on the performance index chosen 
and the type of transfer maneuver considered. Here, we focus 
attention on the minimization of the total characteristic velocity 
AV, which is a measure of the propellant consumption required 
for orbital transfer. 

Problem (IA) . This problem refers to the indirect ascent to 
a 178 NM perigee via a 197 NM apogee. The functional to be minimized 
is given by 

I = AV = AV 22 + AV,, , (7a) 

with 


4V 22 

''a’ /I2r a r ll /(r ll r 22 

+ r 22 n 

- <r a /r 22 )V lC os9 ir 

(7b) 

AV 11 

= V !2r a r 22 /(r U r 22 

+ r^)] 

- V (r a /r ll> • 

(7c) 


Problem (DA) . This problem refers to the direct ascent to 
a 178 NM apogee. The functional to be minimized is given by 

I = AV = AV , (8a) 

with 


AV 


11 


- V ( V r n> - 


,r a /r ll )V l cos1f r 


(8b) 
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Alternative Formulations . Inspection of the performance 
indexes (7) and (8) shows that they include a constant part and 
a part which is linear in the exit value of the horizontal 
component of the inertial velocity = Vcosy. Hence, the 
minimization of AV is equivalent to maximizing the functional 

J = v hl = V 1 cosy 1 . (9a) 

In turn, because of the final conditions (6b) or (6c) , the 
maximization of the functional (9a) implies either maximizing the 
functional 

K = V (9b) 

or minimizing the functional 

H - y 1 . (9c) 

Continuous vs Discrete Control . For both Problems (IA) and 
(DA), two cases are investigated: continuous control and discrete 
control . 

(i) Continuous Control. Here, the bank angle is continuously 
variable. Hence, y (t) is treated as a control variable. 

(ii) Discrete Control. Here, the trajectory is divided into 
s segments along which the bank angle is constant; hence, for 
each segment, the bank angle y^ is treated as a parameter. Because 
the time duration of each segment is also a parameter to be 
optimized, the number of parameters is twice the number of 
segments. The following subcases are studied: s = 2, 3, 4, 5. 

These subcases involve n = 4, 6, 8, 10 parameters respectively . 

ir 
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Algorithm . For both cases (i) and (ii) , Problems (IA) and 
(DA) are optimal control problems of the Bolza type. They can be 
solved using the sequential gradient-restoration algorithm (SGRA) , 
employed in conjunction with the primal formulation (PSGRA, 

Refs . 6-9) . 

The sequential gradient-restoration algorithm is an iterative 
technique which involves a sequence of cycles, each cycle including 
two phases: the gradient phase and the restoration phase. In 
the gradient phase, the value of the augmented functional is 
decreased, while avoiding excessive constraint violation. In 
the restoration phase, the value of the constraint error is 
decreased, while avoiding excessive change in the value of the 
functional. In a complete gradient-restoration cycle, the value 
of the functional is decreased, while the constraints are satisfied 
to a preselected degree of accuracy. Thus, a succession of 
suboptimal solutions is generated, each new solution being an 
improvement over the previous one from the point of view of 
the value of the functional being minimized. 

The convergence conditions are represented by the relations 

P < e x , Q < e 2 - (10) 

Here, P is the norm squared of the error in the constraints? Q 
is the norm squared of the error in the optimality conditions; 
and e,, e 2 are preselected small, positive numbers. 
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5 . Experimental Data 

The following data were used in the numerical experiment 
(for more significant digits, see Ref. 4). 

Earth’s Physical Constants . The radius of the Earth is 
r 0 = 6378 km; the radius ' of the outer edge of the atmosphere is 
r = 6500 km; the thickness of the atmosphere is h = 121.9 km; 

cl 3. 

3 2 

the Earth's gravitational constant is y e = 0.3986E+15 m /sec ; 

the circular velocity at r = r is V = 7.831 km/sec; the angular 

velocity of the Earth is u = 0.7292E-04 rad/sec. 

Spacecraft Data . The mass of the AFE spacecraft is m = 1678 kg; 

2 

the reference surface area is S = 14.31 m ; the spacecraft flies 
at constant angle of attack, a = 17.00 deg; the lift coefficient 
is C T =-0.3707; and the drag coefficient is C = 1.315. 

Atmospheric Model . The assumed atmospheric model is that of 
the US Standard Atmosphere, 1976 (Ref. 20). In this model, the 
values of the density are tabulated at discrete altitudes. For 
intermediate altitudes, the density is computed by assuming an 
exponential fit for the function p (h) . 

Heating Rate . The stagnation point heating rate is assumed 
to be given by the formula 

HR = C/lp/p R ) (V/V r ) 3 - 07 . (ID 

Here, p R is a reference density, the density at h R = 60 km; V R = V a is a 

reference velocity; the constant C represents the stagnation point 

2 

heating rate at p = p and V = V^, and is given by C = 282.3 W/cm ; note 
that the value of this constant is based on a nose radius of one foot. 
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Atmospheric Entry Conditions . In the inertial system, the given 
initial conditions are as follows: the longitude is 0g = -134.52 deg? 
the latitude is = -4.49 deg; the altitude is h Q = 121.9 km; 
the radius is r Q = 6500 km; the velocity is V Q = 10.31 km/sec ; 
the path inclination is y q = -4.49 deg: the heading angle is 

= -28.13 deg; the orbital inclination is i Q = 28.45 deg; the 
longitude of the ascending node is fig = -126.19 deg; the wedge 
angle is rig = 0.00 deg. 

Atmospheric Exit Conditions . In the inertial system, the 
desired final conditions are as follows: the altitude is h^ = 121.9 km; 
the radius is r^ = 6500 km; the orbital inclination is i^ = 28.45 deg; 
the longitude of the ascending node is = -126.19 deg; the 
wedge angle is = 0.00 deg. 

Transfer (IA) . This is the indirect ascent to a 178 NM perigee 
via a 197 NM apogee. The required apogee conditions are as follows: 
the altitude is h 22 = 364.8 km; the radius is r 22 = 6743 km; and 
the path inclination is y 22 =0.00 deg. The required perigee 
conditions are as follows: the altitude is h^ = 329.7 km; the 
radius is r^ = 6708 km; and the path inclination is = 0.00 deg. 

Transfer (DA). This is the direct ascent to a 178 NM apogee. 

The required apogee conditions are as follows: the altitude is 
h lx = 329.7 km; the radius is r ^ = 6708 km; and the path 
inclination is y ^ = 0.00 deg. 
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6 . Numerical Results 

The optimal control problem formulated in Sections 3-4 
was solved with the sequential gradient-restoration algorithm 
for the experimental data outlined in Section 5. Both Transfer 
(IA) and Transfer (DA) were optimized from the point of view of 
the characteristic velocity AV. For each transfer, five optimal 
trajectories (OT) were computed: the continuous control OT; the 
discrete control OT with s = 2 segments; the discrete control OT 
with s = 3 segments; the discrete control OT with s = 4 segments; 
and the discrete control OT with s = 5 segments. For Transfer (IA) , 
a reference trajectory (RT) was also computed: this is a nonoptimal 
discrete-control trajectory with s = 5 segments, obtained by means 
of the restoration algorithm, starting from a nominal trajectory 
supplied in Ref. 5. Thus, a total of 11 trajectories were computed, 
ten optimal trajectories and one reference trajectory. The 
numerical results are shown in Tables 1-2 and Figs. 1-2. 

Table 1 shows the values of the characteristic velocity AV 
and the flight time x for all the trajectories computed. For 
Transfer (IA) , the reference trajectory requires AV = 98.0 m/sec 
and a flight time x = 488 sec; also for Transfer (IA) , the 
optimal trajectories require AV = 81.9 m/sec and a flight time 
x = 796 sec. For Transfer (DA), the optimal trajectories require 
AV = 72.0 m/s ec and a flight time x = 820 sec. To sum up, 
the optimal trajectories are considerably more efficient propulsively 
than the reference trajectory and they are characterized by a 
longer flight time. 
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An interesting result of the analysis is that, for each 
transfer, the optimal continuous-control trajectory and the optimal 
discrete-control trajectories coalesce into a single trajectory: 
a two-subarc trajectory, with the bank angle constant in each 
subarc (bang-bang control) . Specifically, the bank angle 
is near 180 deg in the atmospheric entry phase (positive 
lift projection phase) and is near 0 deg in the atmospheric 
exit phase (negative lift projection phase) . This is true for 
both Transfer (IA) and Transfer (DA). See Figs. 1-2. 

In the light of the above statement, we now focus attention 
on only three trajectories: for Transfer (IA) , the reference 
trajectory and the two-subarc optimal trajectory; for Transfer 
(DA), the two-subarc optimal trajectory. For these trajectories, 
Table 2 shows the values of the following quantities: the 
characteristic velocity; the flight time; the minimum altitude; 
the peak dynamic pressure; the peak heating rate; the peak change 
of the orbital inclination; the peak change of the longitude of 
the ascending node; and the peak value of the wedge angle. 

Table 2 illustrates clearly the advantages of the OTs 
with respect to the reference trajectory. Concerning the longitudinal 
motion, by comparison with the RT, the OTs are characterized by 
smaller AV and also by smaller values of the peak dynamic pressure 
and the peak heating rate; this is due to the fact that the 
minimum altitude of the OTs is about 3 km higher than the minimum 
altitude of the RT. Concerning the lateral motion, by 
comparison with the RT, the OTs exhibit smaller values of the 
peak change of the orbital inclination and the peak change 
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of the longitude of the ascending node; hence, the peak wedge 
angle of the OTs is smaller than the corresponding quantity for 
the RT. Indeed, it is suprising that max|Ai| and max | Aft | are 
nearly zero for the OTs;hence, max(n) is nearly zero for the OTs. 
This means that, for efficient flight, the motion of the spacecraft 
must be nearly planar in an inertial space; in other words, one 
must avoid energy dissipation associated with the lateral motion. 

More details on the trajectories computed can be found in 
Figs. 1 and 2. Figure 1 refers to Transfer (IA),and Fig. 2 refers 
to Transfer (DA) . These figures include the time histories of 
the following quantities: the bank angle (Figs. 1A and 2A) ; the 
altitude (Figs. IB and 2B) ; the relative velocity (Figs. 1C and 2C) ; 
the relative path inclination (Figs. ID and 2D); the dynamic 
pressure (Figs. IE and 2E) ; the heating rate (Figs. IF and 2F) ; 
and the wedge angle (Figs. 1G and 2G) . 

Remark . With reference to the two-subarc OTs of Transfer (IA) 
and Transfer (DA) , we observe that the value of the bank angle 
is y = 176.7 deg for the atmospheric entry phase and y = 5.5 deg 
for the atmospheric exit phase. Let L denote the lift, and let 

L = Lcosy, = Lsiny (12) 

P <3 

denote the vertical component and the lateral component of the 

lift. Let C T , C T , C denote the associated lift coefficients, 

Li Lp LQ 


and observe that 
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C T = CLcosii, C T = CLsiny. (13) 

J_iP Jj J-iCJ -Li 

We recall that the AFE spacecraft flies at constant angle of 

attack a = 17 deg, corresponding to CL = -0.3707 and CL = 1.315. 

-Li L) 

With this understanding, we see that, for the atmospheric entry 
phase , 

C T = -0.998 CL , C = 0.058CL; (14) 

Lp L Lq L 

and we see that, for the atmospheric exit phase, 

CL = 0.995 C T , CL = 0.096 C T . (15) 

Lp L' Lq L 

Therefore, the following comments arise: 

(i) in the atmospheric entry phase, the vertical component 
of the lift is directed upward and its modulus is only 2/1000 
smaller than the lift modulus; this causes the path inclination 
to increase gradually from the entry negative value to nearly 
zero value; 

(ii) in the atmospheric exit phase, the vertical component 
of the lift is directed downward and its modulus is only 5/1000 
smaller than the lift modulus; this offsets the centrifugal force 
effects due to the curvature of the Earth, so as to ensure exit 
conditions compatible with the desired apogee requirement; 

(iii) the lateral component of the lift during the atmospheric 
entry phase and the lateral component of the lift during the 
atmospheric exit phase have the same sign and the same order of 
magnitude; their modulus is 5.8% of the lift modulus in the 
atmospheric entry phase and 9.6% of the lift modulus in the atmospheric 
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exit phase; these lateral lift components are directed in such a 
way that they nearly offset the effects due to the Earth's 
rotation, so as to ensure exit conditions compatible with the 
desired wedge angle requirement. 
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7 . Guidance Implications 

The results of Section 6 provide a foundation for developing 
a more efficient guidance scheme for the AFE spacecraft. For 
information concerning the existing guidance scheme, see Refs. 5 
and 21. 


First, we recall the basic idea of the existing guidance 
scheme. For the AFE spacecraft, the angle of attack varies from 
7 to 27 deg, and the lift coefficient varies from -0.47 to -0.21. 
Therefore, if the AFE spacecraft were to be controlled via only 
the angle of attack, the lift coefficient range would be | AC^ | = 0.26 
To offset the above difficulty, the existing guidance scheme 
employs a constant angle of attack a = 17 deg, corresponding to 
C T = -0.37, coupled with a variable angle of bank. 

Let Lp = Lcosy denote the projection of the lift vector in 
the vertical direction; let C T = C T cosy denote the projected 

LiP J-i 

lift coefficient. Because C T is constant and cosy varies from 

1j 

-1 to +1, the projected lift coefficient range is | AC | =0.74. 

This is nearly triple the lift coefficient range |AC l | = 0.26. 

To sum up, for the AFE spacecraft, the existing guidance 
scheme employs constant angle of attack and variable angle of 
bank in order to achieve the following goals: for the longitudinal 
motion, to dissipate excess velocity so that the spacecraft, 
upon exiting the atmosphere, can ascend to a specified apogee; 
and for the lateral motion, to keep the instantaneous orbital plane 
close to the initial orbital plane, hence to keep the wedge 
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angle close to zero. These goals are achieved as follows: for 
the longitudinal motion, one controls the modulus of the bank 
angle; for the lateral motion, one controls the sign of the bank 
angle; thus, whenever the wedge angle exceeds a specified bound, 
the sign of the bank angle is reversed. 

Unquestionably, the existing guidance scheme has obvious 
merits: (a) simplicity; and (b) the fact that the projected lift 

coefficient range is three times the lift coefficient range. 

However, there are weaknesses associated with the control of the 
lateral motion: 

(A) There is no equilibrium point for the control of the 
lateral motion. When n = 0, it would be desirable to have n = 0, 
so that the wedge angle continues to be zero. However, y is 
determined from the requirement of longitudinal motion control; 
hence, generally speaking, n ^ 0 when ri = 0 . This means that 
the system is not stable around y = 0. 

(B) There is a bank angle error due to noninstantaneous 
switches. The implementation of the existing guidance scheme ideally 
requires instantaneous switches in sign(y), so as to keep cosy 
unchanged. In practice, this is not possible, since |y| <^20 deg/sec. 
For example, a switch from y = + 170 to y = -170 deg requires 

At = 17 sec. If the correct y-values are y = +170 and y = -170 deg, 
this means that, during the time interval At = 17 sec, there is 
an error in the value of y which is required for the control of the 
longitudinal motion. 
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(C) There is a contradiction between accuracy and stability. 

For accurate control of the lateral motion, the wedge angle tolerance 

should be small. On the other hand, if the wedge angle tolerance 

is small, the number of switches in sign(y) increases; therefore, 

the system stability becomes worse. 

To offset the above difficulties, consideration should be 

given to developing a new guidance scheme, based on the properties 

of the optimal trajectories. The basic ideas of the new guidance 

scheme are: to control both the angle of attack (near 17 deg) and 

the angle of bank; and to utilize a single switch of the bank 

angle, consistently with the optimal trajectory properties, so 

as to have L > 0 in the entry phase and L < 0 in the exit phase; 

P P 

here, L is the projected lift. The new guidance scheme retains 
P 

the basic advantage of the existing guidance scheme, in that 
the projected lift coefficient range is three times the lift 
coefficient range. In addition, it has the following advantages: 

(A) There is an equilibrium point for the control of the 
lateral motion, due to the fact that both the angle of attack and 
the angle of bank are being varied. Hence, when n = 0, it 

is possible to have ri = 0, which means that the system is stable 
around n = 0. 

(B) There is less bank angle error due to noninstantaneous 
switches. For the existing guidance scheme, the effects due to 
the noninstantaneous switch from +y to — y can be serious, since 
there are multiple switches and they occur randomly. For the new 
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guidance scheme, the effects due to the noninstantaneous switch 
from +y to -y are less serious, owing to the fact that there 
is only one switch, which occurs whenever a preselected velocity 
is achieved. 

(C) There is improved accuracy and stability in the control 
of the lateral motion. This is because the lateral motion is no 
longer controlled by the sign of the bank angle, but by the 
values of both the bank angle and the angle of attack. 

(D) The new guidance scheme is based on the properties of 
the optimal trajectories. Hence, it preserves the good properties 
of the optimal trajectories concerning the characteristic velocity, 
the peak dynamic pressure, the peak heating rate, and the peak 
wedge angle. 

(E) The new guidance scheme uses the same hardware (AFE 
configuration, measurements, sensors, and reaction control system) 
as the existing guidance scheme. The only difference is in the 
software (computer code). Therefore, if acceptable, a change from 
the existing guidance scheme to the new guidance scheme would be 


of limited cost. 
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8 . Conclusions 

In this report, the trajectories of an AFE spacecraft are 
analyzed in a 3D-space, employing the full system of 6 ODEs 
describing the atmospheric pass. It is assumed that the angle 
of attack is constant and that the spacecraft is controlled via 
only the angle of bank. It is also assumed that the atmospheric 
eptry conditions are given and that the atmospheric exit conditions 
are adjusted in such a way that the following conditions are 
satisfied: (a) the atmospheric velocity depletion is such that, 
after exiting, the AFE spacecraft first ascends to a specified 
apogee and then descends to a specified perigee; and (b) the 
exit orbital plane is identical with the entry orbital plane. 

Under the above conditions, optimal trajectories are 
determined by minimizing the total characteristic velocity, hence 
the propellant consumption required for orbital transfer. The 
resulting optimization problems are solved by means of the sequential 
gradient-restoration algorithm. 

Two possible transfers are considered: (IA) indirect ascent 

to a 178 NM perigee via a 197 NM apogee; and (DA) direct ascent 
to a 178 NM apogee. For both transfers, two cases are investigated: 
(i) the bank angle is continuously variable; and (ii) the trajectory 
is divided into segments along which the bank angle is constant. 

For case (ii) , the following subcases are studied: two segments, 
three segments, four segments, and five segments; because the time 
duration of each segment is optimized, the above subcases involve 
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four, six, eight, and ten parameters, respectively. For comparison 
purposes and only for Transfer (IA) , a reference trajectory is 
also considered: this is a five-segment trajectory, close to the 
nominal trajectory given in Ref. 5. 

From the extensive numerical computations, the following 
conclusions arise: 

(i) The optimal trajectories are two-subarc trajectories, 
with the bank angle constant in each subarc; hence, the control 
is bang-bang. 

(ii) In the atmospheric entry phase, the bank angle is 
near 180 deg, yielding a positive projected lift, which in turn 
causes the path inclination to increase gradually from the entry 
negative value to nearly zero value. 

(iii) In the atmospheric exit phase, the bank angle is near 
0 deg, yielding a negative projected lift, which offsets the 
centrifugal force effects due to the curvature of the Earth, so 
as to ensure exit conditions compatible with the desired apogee 
requirement. 

(iv) The lateral component of the lift during the atmospheric 
entry phase and the lateral component of the lift during the 
atmospheric exit phase have the same sign and the same order of 
magnitude; they are directed in such a way that they nearly offset 
the effects due to the Earth's rotation. In this way, the 
instantaneous orbital plane is almost identical with the initial 
orbital plane, meaning that the wedge angle is nearly zero during 
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the atmospheric pass. This means that, for efficient flight, the 
motion of the AFE spacecraft is nearly planar in an inertial space; 
in other words, one must avoid energy dissipation associated with 
the lateral motion. 

(v) Comparison of the optimal trajectories and the reference 
trajectory shows that the OTs are superior to the RT in terms 

of the main quantities of interest, namely, the characteristic 
velocity, the peak dynamic pressure, the peak heating rate, and 
the peak wedge angle. In particular, for Transfer (IA) , the 
characteristic velocity is AV = 98.0 m/sec for the RT and 
AV = 81.9 m/sec for the OT. For Transfer (DA) , the characteristic 
velocity of the OT is 72.0 m/sec . 

(vi) The results on optimal trajectories have important 

guidance implications. They suggest the idea of developing a 

new guidance scheme based on the following precepts: to control 

both the angle of attack (near 17 deg) and the angle of bank; 

and to utilize a single switch of the bank angle, so as to have 

L > 0 in the entry phase and L < 0 in the exit phase. This 

p P 

new guidance system is promising, and it appears to be an improvement of 
the existing guidance scheme in the control of the lateral motion 
of the AFE spacecraft. 
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Table 1. Characteristic velocity and flight time. 


Trajectory 

Control 

s 

AV (m/sec) 

t (sec) 

IART 


Discrete 

5 

97.97 

487.6 

IAOT 


Continuous 

00 

81.86 

795.3 

IAOT 


Discrete 

2 

81.86 

795.7 

IAOT 


Discrete 

3 

81.86 

795.7 

IAOT 


Discrete 

4 

81.86 

795.8 

IAOT 


Discrete 

5 

81.84 

796.2 

DAOT 


Continuous 

CO 

72.05 

817.6 

DAOT 


Discrete 

2 

72.03 

819.5 

DAOT 


Discrete 

3 

72.03 

819.6 

DAOT 


Discrete 

4 

72.03 

819.6 

DAOT 


Discrete 

5 

72.01 

819.9 

IA = 

indirect ascent, DA = 

direct ascent. 


OT = 

optimal 

trajectory, RT = 

reference 

trajectory. 


S = 

number 

of segments. 
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Table 2. Comparison of AFE trajectories. 


Quantity 

IART 
s = 5 

IAOT 
s = 2 

DAOT 
s = 2 

Units 

AV 

97.97 

81.86 

72.03 

m/sec 

T 

487.6 

795.7 

819.5 

sec 

min (h) 

74.6 

77.5 

77.5 

km 

max (DP) 

1629 

1174 

1174 

N/m 2 

max (HR) 

158 

148 

148 

W/cm 2 

max | Ai | 

1.08 

0.01 

0.01 

deg 

max | Aft | 

0.32 

0.00 

0.00 

deg 

max (n) 

1.08 

0.01 

0.01 

deg 

IA = indirect 

ascent , 

DA = direct 

ascent , 



OT = optimal trajectory, RT = reference trajectory, 
s = number of segments. 
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Abstract . Tft«-co pr>rt. in t -hia t hi xd a . s^± & s (p r fc. l~4t 
dg^lirfig w-ith the dete rm i nat i on of op t i mal — fra jpct -o r i e e— ftn — the a e r e- 

H rtffrV 

assisted flight experiment (AFE^. The intent of this experiment is to 
simulate a GEO— to-LEO transfer , where GEO denotes a Geosynchronous 
Earth orbit and LEO denotes a low Earth orbit. Specif ica - ll y , 

^-fche— AFE s pac e craft is r -&Xeased— £rom-^-th&- -spaee- shuttle -arrd 


-fee— a c hiev e 


— e ntry condi tions iden Lical — wit- h- t - hc se — cf 


A- s p acecra ft r eturning fr o m G EO. During-tiie ePmOspher ic~ pass , 

-fe iie cU ' i gle of aTFaCk— i-9— k e pt -eon s tant, — and— the -angle of bank- ie- 
eontrcdrled-in- stroh~a-way— that ^the following conditions are satisfied: 
(a) the atmospheric velocity depletion is such that, after 
exiting, the AFE spacecraft first ascends to a specified apogee and 
then descends to a specified perigee; and (b) the exit orbital 
plane is identical with the entr^orbital plane The final maneuver ; 
not analyzed here, includes the rendezvous with and the capture 
by the space shuttle. 

In. this r e port ^ the trajectories of an AFE spacecraft are 
analyzed in a 3D-space, employing the full system of 6 ODEs describing 
the atmospheric pass. The atmospheric entry conditions are given, 
and the atmospheric exit conditions are adjusted in such a way that 


ac qu i rement s- 




the total characteristic velocity, hence the propellant consumption 

* A 

£.r<_ ** i n * ** < 2 . 

required for orbital transfer^ Two possible transfers are 
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considered: (IA) indirect ascent to a 178 NM perigee via a 197 NM 

apogee; and (DA) direct ascent to a 178 NM apogee. 

For both transfers, two cases are investigated: (L) the 

bank angle is continuously variable; and (Q- ) the trajectory is 

divided into segments along which the bank angle is constant. 

For case (:$.), the following subcases are studied: two segments, 

three segments, four segments, and five segments ; because the time 

duration of each segment is optimized, the above subcases involve 

four, six, eight, and ten parameters, respectively. 

•x«. aUc. cuvA.lijV’te om.I'Js.’om , 

•ft — s urprising result of tho analysis is fe-ha ^ the optimal 
trajectories of cases (i) and (A) a coalesce into a single trajectory: 
a two-subarc trajectory, with the bank angle constant in each 
subarc -^b ang - ba - ny -c oivtrul j-. Specifically, the bank angle is 
near 180 deg in the atmospheric entry phase 'tpositive-lift prc>3 ei^te-ion- 
_pka*reT and is near 0 deg in the atmospheric exit phase 

Li v e lift piojection phnsc T'l -Another sur pri s ing r e sult is that > 
during the atmospheric pass, the peak values of the changes of the 
orbital inclination and the longitude of the ascending node are 
nearly zero; hence, the peak value of the wedge angle (angle between 
the instantaneous orbital plane and the initial orbital plane) is 
nearly zero. This means that the motion of the spacecraft is 
nearly planar in an inertial space. 

The guidance implications of the above results are discussed. 

Key Words. Flight mechanics, hypervelocity flight, atmospheric 
flight, optimal trajectories, aeroassisted flight experiment, 
aeroassisted orbital transfer, guidance. 


